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We introduce a classification scheme of the post-merger dynamics and gravitational-wave emission 
in binary neutron star mergers, after identifying a new mechanism by which a secondary peak in the 
gravitational-wave spectrum is produced. It is caused by a spiral deformation, the pattern of which 
rotates slower with respect to the double-core structure in the center of the remnant. This secondary 
peak is typically well separated in frequency from the secondary peak produced by a nonlinear 
interaction between a quadrupole and a quasi-radial oscillation. The new mechanism allows for 
an explanation of low-frequency modulations seen in a number of physical characteristics of the 
remnant, such as the central lapse function, the maximum density and the separation between the 
two cores, but also in the gravitational-wave amplitude. We find empirical relations for both types 
of secondary peaks between their gravitational-wave frequency and the compactness of nonrotating 
individual neutron stars, that exist for fixed total binary masses. These findings are derived for equal- 
mass binaries without intrinsic neutron-star spin analyzing hydrodynamical simulations without 
magnetic field effects. Our classification scheme may form the basis for the construction of detailed 
gravitational-wave templates of the post-merger phase. We find that the quasi-radial oscillation 
frequency of the remnant decreases with the total binary mass. For a given merger event our 
classification scheme may allow to determine the proximity of the measured total binary mass to 
the threshold mass for prompt black-hole formation, which can, in turn, yield an estimate of the 
maximum neutron-star mass. 

PACS numbers: 04.30.Tv,26.60.Kp,97.60.Jd,04.40.Dg 


I. INTRODUCTION 

Neutron star (NS) mergers are strong emitters of grav¬ 
itational waves (GWs) and thus among the prime targets 
for the upcoming GW detectors Advanced LIGO [l| and 
Advanced Virgo [^. Future GW observations of such 
events Q could reveal the properties of high-density mat¬ 
ter and NSs (see e.g. [J-Q for reviews). The merger 
will likely result in the formation of a differentially 
rotating, strongly oscillating remnant [ol-l^. Specifi¬ 
cally, it has been shown that detections of its domi¬ 
nant post-merger GW frequency /peak would strongly 
constrain the radius and the maximum mass of nonro¬ 
tating NSs [H, [ 2 ^, m, m, . Additionally, there exist 
potent ially detectable secondar y pe aks at lower frequen¬ 
cies [i,[Il[Il,[ll[Il[ll,[2i,[2l[2BM[3l[3i-l4ll: which 
could improve the constraints on NS properties if the 
mechanism by which they appear is understood. 

In this study we find that two different mechanisms 
are at work for producing the low-frequency secondary 
peaks in the GW spectrum at frequencies below the main 
peak of the post-merger phase. Specifically, we show that 
apart from a nonlinear combination frequency , there 
exists also a distinct mechanism that generates a sec¬ 
ondary GW peak by the rotating pattern of a deformation 
of spiral shape. This deformation is initially produced at 
the time of merging and is then sustained for a few ro¬ 
tational periods. The consideration of the two different 
mechanisms leads to a unified classification scheme for 
the post-merger dynamics and GW emission. For high- 
mass binaries (relative to the threshold mass to prompt 


black-hole collapse), the nonlinear combination frequency 
dominates, while for low-mass binaries it is the spiral de¬ 
formation that produces the strongest secondary peak. 
Both are simultaneously present and can produce peaks 
of comparable strength for intermediate binary masses. 
Hence, the secondary peaks cannot simply be explained 
by the dynamics of the double-core structure alone, which 
forms the inner part of the remnant, as in [36l - l3^ . Iden¬ 
tifying the importance of the rotating deformation, we 
explain for the first time the existence of a low-frequency 
modulation appearing in addition to the quasi-radial os¬ 
cillation mode in several physical quantities characteriz¬ 
ing the remnant. 

After clarifying the nature of the secondary peaks in 
the GW spectrum, we find mass-dependent empirical re¬ 
lations for the different peaks as a function of NS com¬ 
pactness. We rule out the existence of a universal, mass- 
independent and EoS independent relation for the sec¬ 
ondary peaks, a relation recently proposed in [Hill. 
Such a relation can only be found within a small sample 
of EoSs for a very limited mass range (not the same for 
all EoSs). Our new relations can be used to optimize the 
constraints on the EoS in the case of the simultaneous 
detection of several post-merger frequencies. The under¬ 
standing of the most prominent features of the GW spec¬ 
trum is a prerequisite for constructing GW templates for 
the post-merger phase, which could enhance the detec¬ 
tion prospects compared to unmodelled searches [13, |4l| 
for the Advanced LIGO and Advanced Virgo detectors 
and their discussed upg rades [4^ - l43| . Eor the planned 
Einstein Telescope direct detections of secondary 
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peaks are a viable prospect [H, [13,113, IHl • 


II. NATURE OF SECONDARY GW PEAKS 

We investigate mergers of equal-mass, intrinsically 
non-spinning NSs with a 3D relativistic smoothed par¬ 
ticle hydrodynamics (SPH) code, which imposes the con¬ 
formal flatness condition on the spatial metric [13, [13 
to solve Einstein’s field equations and incorporates en¬ 
ergy and an gula r momentum losses by a GW backreac- 
tion scheme |l3 . [11] (see [H, [H, [H, [13, [13] for details on 
the code, the setup, resolution tests and model uncertain¬ 
ties). Comparisons to other numerical setups and also 
models with an approximate consideration of neutrino ef¬ 
fects show an agreement in determining the post-merger 
spectrum within a few per cent in the peak frequen¬ 
cies [23 - E3 . [H, [H-®!- Magnetic field effects are neg¬ 
ligible for not too high initial field strengths [11| . We ex¬ 
plore a representative sample of ten microphysical, fully 
temperature-dependent equations of state (EoSs) (see 
Table I in (sl] and Fig. [5] in this work for the mass-radius 
relations of non-rotating NSs of these EoSs) and consider 
total binary masses Mtot between 2.4 Mq and 3.0 Mq. 
In this work we consider only NSs with an initially ir- 
rotational velocity profile because known spin periods in 
observed NS binaries are slow compared to their orbital 
motion (see e.g. 0) , and simulations with initial intrin¬ 
sic NS spin suggest an impact on the post-mer ger features 
of the GW signal only for very fast spins [l3 [35l. [^ . 

First, we focus on a reference model for the moderately 
stiff DD2 EoS [U, [53| with an intermediate binary mass 
of Mtot = 2.7 Mq. Figured] shows the x-polarization of 
the effective amplitude /leff.x = ^x(/) • / (with being 
the Fourier transform of the waveform /ix) vs. frequency 
/ (reference model in black). Besides the dominant /peak 
frequency [^ . there are two secondary peaks at lower 
frequencies (/ 2-0 and /spiral) with comparable signal-to- 
noise ratio. Both are generated in the post-merger phase, 
which can be seen by choosing a time window covering 
only the post-merger phase for computing the GW spec¬ 
trum. 

The secondary peak shown as / 2-0 is a nonlinear com¬ 
bination frequency between the dominant quadrupolar 
/peak oscillation and the quasi-radial oscillation of the 
remnant, as described in |25l |. We confirm this by per¬ 
forming additional simulations, after adding a quasi- 
radial density perturbation to the remnant at late times. 
The frequency /o of the strongly excited quasi-radial os¬ 
cillation is determined by a Fourier analysis of the time- 
evolution of the density or central lapse function and co¬ 
incides with the frequency difference /peak — / 2 -o- As 
in [^, the extracted eigenfunction at /o confirms the 
quasi-radial nature. 

The secondary /spiral peak is produced by a strong de¬ 
formation initiated at the time of merging, the pattern 
of which then rotates (in the inertial frame) slower than 
the inner remnant and lasts for a few rotational peri- 
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FIG. 1: GW spectra of 1.35-1.35 Mq mergers with the 
DD2 [si],[53 (black), NL3 [53,[53| (blue) and LS220 [H (red) 
EoS (cross polarization along the polar axis at a reference dis¬ 
tance of 20 Mpc). Dashed lines show the anticipated unity 
SNR sensitivity curves of Advanced LIGO [3 (red) and of the 
Einstein Telescope [43 (black). 

ods, while diminishing in amplitude. Figured shows the 
density evolution in the equatorial plane, in which one 
can clearly identify the two antipodal bulges of the spi¬ 
ral pattern, which rotate slower than the central parts 
of the remnant. In this early phase the inner remnant is 
still composed of two dense cores rotating around each 
other (this is the nonlinear generalization of an m = 2 
quadrupole oscillation producing the dominant /peak)- 
Extracting the rotational motion of the antipodal bulges 
in our simulations, we indeed find that their frequency 
equals /spirai/2 producing gravitational waves at /spiral 
(compare the times in the right panels in Fig. [2] recall 
the factor two in the frequency of the GW signal com¬ 
pared to the orbital frequency of orbiting point particles). 
In Fig. 12] the antipodal bulges are illustrated by selected 
fluid elements (tracers), which are shown as black and 
white dots, while the positions of the individual centers 
of the double cores are marked by a cross and a circle. 
(We define the centers of mass of the double cores by 
computing the centers of mass of the innermost 1000 
SPH particles of the respective initial NSs and then fol¬ 
lowing their time evolution.) While in the right panels 
the antipodal bulges completed approximately one orbit 
within one millisecond (« -r-^—), the double cores moved 

' /spiral ' 

further ahead, i.e. with a significantly higher orbital fre¬ 
quency. Examining the GW spectrum and considering 
different time intervals, we find that the presence of the 
/spiral peak agrees with the appearance and duration of 
the spiral deformation of the remnant. 

In the upper right panel of Fig. 2, the spiral deforma¬ 
tion can be seen to initially reach deep inside the rem¬ 
nant. We approximately determine the amount of matter 
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FIG. 2: Rest-mass density evolution in the equatorial plane for the 1.35-1.35 Mq merger with the DD2 EoS (rotation counter¬ 
clockwise). (The rest-mass density is shown with a variable linear scale relative to pmax. A low number of contour levels 
is chosen for illustrative reasons; the underlying simulation data is smoother than it appears with the chosen color coding.) 
Black and white dots trace the positions of selected fluid elements of the antipodal bulges, which within approximately one 
millisecond complete one orbit (compare times of the right panels). The orbital motion of this pattern of spiral deformation 
produces the /spUai peak in the GW spectrum at 2/(1 ms) (Fig. 1). The cross and the circle mark the double cores, which 
rotate significanty faster than the antipodal bulges represented by the dots (compare times of the different panels). 


which belongs to the two antipodal bulges that are ro¬ 
tating slower compared to the double cores. This matter 
amounts to several tenths of Mq and is thus sufficient 
to explain the strength of the /spiral GW peak. In ad¬ 
dition, we find that the /spiral GW peak can be roughly 
reproduced in a toy model, where the two bulges orbit as 
point particles around the central double-core structure 


for a duration of a few milliseconds. Note that this toy 
model differs significantly from the one in [3^ , which con¬ 
siders only the two cores to be contributing to the GW 
signal and considers only a single instantaneous orbital 
frequency of the system. 

Furthermore, we take advantage of the quadrupole for¬ 
malism to compute GW spectra considering only certain 
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FIG. 3: Early post-merger spectrum of 1.35-1.35 Mq merger 
with the DD2 EoS (black) and GW spectra computed for fluid 
elements with densities above or below 50% of the maximum 
density. 

parts of the remnant, which are defined by using either a 
density or a spatial cut-off. In Fig.jSjwe demonstrate that 
the dominant quadrupole fpeak frequency is generated 
mainly by those regions of the remnant which encom¬ 
pass densities exceeding 50% of the instantaneous maxi¬ 
mum density Pmax- In contrast, most power of the /spiral 
peak originates from densities below O.Sprnax, which cor¬ 
responds to the outer parts of the remnant, where the two 
bulges form (see Fig.|2|). Similar conclusions are reached 
when a spatial cut-off instead of a density cut-off is used. 

In models where /spiral dominates over / 2-0 (see our 
classification below), the presence of the two rotating an¬ 
tipodal bulges explains the appearance of a particular 
low-frequency modulation with /peak — /spiral seen in the 
time evolution of the central lapse function (blue curve 
in Fig. 2]), of pmax, of the size of the central remnant and 
of the separation between the two cores. The same mod¬ 
ulation occurs as a beat frequency in the time evolution 
of the GW amplitude (see e.g. Fig. 1 in El)- The low- 
frequency modulation with /peak ~ /spiral can coexist with 
the quasi-radial oscillation, which has a higher frequency 
/o. 

The above modulation frequency is associated with the 
rotating bulges and is explained as follows: the central 
remnant forms an elongated structure, around which the 
two bulges rotate, lagging behind. The characteristics 
of the remnant are modulated depending on the orien¬ 
tation of the antipodal bulges with respect to the dou¬ 
ble cores: the compactness is smaller, the central lapse 
function larger and the GW amplitude maximal when 
the bulges and the cores are aligned (lower right panel 
in Fig. HI). Instead, when the bulges and the cores are 
orthogonal to each other (lower left panel), the compact¬ 
ness is largest, the central lapse function smaller and the 



FIG. 4: Evolution of the central lapse function for 1.35- 
1.35 Mq mergers with the NL3 (blue) and LS220 (red) EoS 
(time-shifted to the same merging time). The evolution of 
models with an added quasi-radial perturbation is shown 
(green). For Type II and Type III mergers, like the model 
with the NL3 EoS shown here, one recognizes a low-frequency 
oscillation in addition to the quasi-radial mode, which is also 
strongly excited in the perturbed model at late times. 


GW amplitude small. 

In a frame corotating with the central remnant (i.e. 
a frame rotating with frequency ~ /peak/2 w.r.t. the 
inertial frame) the two bulges counter-rotate with a fre¬ 
quency of A/ = (/peak - /spirai)/2 which corresponds to 
the /peak — /spiral modulation in various quantities (no¬ 
tice the TT-symmetry of the system, which compensates 
the factor 1/2). 

Our simulations also show that the rotational fre¬ 
quency of the centers of the double-core structure, al¬ 
though not constant, significantly exceeds /spirai/2 at any 
time. Hence, variations in the angular frequency of the 
double-core structure alone are not sufficient to inter¬ 
pret the secondary peaks, and the simultaneous presence 
of both /peak and /spiral cannot simply be attributed to 
a single instantaneous angular frequency of the system, 
as suggested in [s^-js^. We note that the simultane¬ 
ous presence of two (orbital) frequencies (of the double¬ 
core structure and of the antipodal bulges) naturally ex¬ 
plains strong time variations of the instantaneous GW 
frequency as seen in many simulations, e.g. (^. Is^. 

Also, a single, initially strongly varying instantaneous 
frequency as explanation for the peaks in the spectrum 
is incompatible with the fact that the most pronounced 
peak in the spectrum occurs already at early times. 
Initially, the instantaneous frequency strongly oscillates 
around the frequency of the dominant peak. According 
to the picture of , peaks form at frequencies at which 
the instantaneous frequency spends most time, i.e. at the 
extrema of the instantaneous frequency. In this case, one 
would thus not expect that the dominant /peak is strong 
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at early times. In contrast, in our simulations the domi¬ 
nant /peak has a substantial strength if one considers only 
the first few milliseconds. 


III. CLASSIFICATION OF POST-MERGER 
DYNAMICS AND GW EMISSION 

We have applied the above analysis tools (GW spectra, 
determination of /o from perturbed models, rotational 
frequency of the antipodal bulges and of the double cores, 
GW spectra with different cut-off densities) for a num¬ 
ber of representative models, varying the binary mass 
and stiffness of the EoS. The results fully confirm the 
generic picture described above. Considering this larger 
set of models, we identify, based on the relative strength 
between / 2-0 and /spiral, three different types of post¬ 
merger dynamics and GW spectra for remnants which 
survive for more than several milliseconds. 

• Type I: When the total binary mass Mjot is not too 
far from the threshold mass for prompt quasi-radial 
collapse of the remnant for a given EoS [32|, the 
evolution of the central lapse function (and of Pmax) 
is dominated by a very strong quasi-radial oscilla¬ 
tion of the remnant, see lower curve in Eig. 21 For 
such models the two initial NSs are more centrally 
condensed and they merge with higher impact ve¬ 
locity (Fig. 3 in [53). Because of the strongly ex¬ 
cited quasi-radial oscillation, / 2-0 is the strongest 
secondary peak in the GW spectrum, while /spiral 
is much weaker, likely because for more compact 
NSs the formation of the spiral pattern is less pro¬ 
nounced. There can be a partial overlap between 
/ 2 _o and /spiral, see red curve in Fig. [TJ 

• Type II: For intermediate total binary masses, / 2-0 
and /spiral have a comparable strength in the GW 
spectrum and the two types of secondary peaks are 
well separated (see black curve Fig. 1). This is 
the generic type, which we discussed earlier in our 
reference model. In characteristic quantities such 
as the central lapse function the quasi-radial os¬ 
cillation frequency /o as well as the low-frequency 
modulation with /peak — /spiral are clearly noticable. 

• Type III: When the total binary mass Mtot is sig¬ 
nificantly below the threshold mass for quasi-radial 
collapse, the time evolution of the central lapse 
function (as well as of Pmax and of the radius of the 
remnant) is dominated by the /peak — /spiral modu¬ 
lation that we explained in the previous section as a 
result of the rotating spiral pattern with the two an¬ 
tipodal bulges. In the evolution of the central lapse 
function this modulation typically has a smaller 
amplitude than Type I variations (Fig. |4]). The 
quasi-radial oscillation /o is also present, but with 
much smaller amplitude than the dominant modu¬ 
lation because the smaller NS compactness implies 



FIG. 5: Different types of post-merger dynamics and GW 
emission of the different merger models visualized by the 
mass-radius relations of non-rotating NSs of the EoSs con¬ 
sidered in this work. The outcome of a given calculation with 
Mtot is shown as symbol at Mtot/2 plotted on the mass-radius 
relation of the EoS employed in the simulation. Red squares 
indicate Type I, black crosses stand for Type II, and blue cirl- 
ces mark Type III. See text for definitions of different types 
of post-merger dynamics and GW emission. 

a smaller impact velocity and thus a weaker exci¬ 
tation of the quasi-radial oscillation. The smaller 
NS compactness also allows for a stronger spiral de¬ 
formation. Consequently, the dominant secondary 
peak in the GW spectrum is /spiral, while / 2-0 is 
either very weak or hidden inside the background 
(see blue curve in Fig. [T|). 

For a given EoS there is a continuous transition be¬ 
tween the different types of post-merger dynamics de¬ 
pending on the total binary mass. Types I and III 
are the limiting cases of the more generic Type II. No¬ 
tice that since the threshold for quasi-radial collapse is 
EoS-dependent [s^, the different types cover a differ¬ 
ent mass range for each EoS. For a total binary mass 
of Mtot = 2.7 Mq all three types are possible depending 
on the EoS, where very soft EoSs yield Type I mergers 
and very stiff EoSs lead to Type III dynamics (see GW 
spetra in Fig. [1}. This is also shown in Fig. [SJ which 
provides an overview of the types of the different models 
considered in this study. The type of a given simulation 
is indicated by a symbol plotted at the mass of the in¬ 
dividual inspiralling stars on the mass-radius relation of 
the EoS which was used in the calculation (e.g. the re¬ 
sults of the 1.2-1.2 Mq binaries are displayed at the radii 
of NSs with 1. 2 Mq). As described, there is a continu¬ 
ous transition between the different classes, which is why 
one should consider the transitions between the different 
types in Fig. [5] as being tentative and a slightly differ¬ 
ent classification may be possible at the borders between 
the different types. More quantitative definitions of the 
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FIG. 6: /peak, /spiral and / 2-0 for mergers with ten different 
EoSs and Mtot=2.7 M© vs. the compactness M/R for nonro¬ 
tating, single NSs. Solid lines show empirical relations. The 
dashed line is taken from (see text for explanations). 


types may be useful in the future. Still, one can clearly 
identify a diagonal band of Type II mergers for interme¬ 
diate binary masses, and also the binary setups leading 
to the limiting cases of Type I or Type III are seen to 
form roughly diagonal bands. 

For 2.4 Mq < Mtot < 3.0 Mq we find that /spiral 
typically ranges between /peak - 0.5 kHz and /peak - 
0.9 kHz, while / 2-0 ranges between /peak — 0.9 kHz and 
/peak — 1.3 kHz. This property will be useful for iden¬ 
tifying either / 2-0 or /spiral (or both) in future GW ob¬ 
servations. Furthermore, we find that /peak — / 2 -o(= /o) 
decreases with increasing Mtot in all models for which 
/ 2-0 is clearly present, in agreement with the fact that 
the quasi-radial frequency decreases near the threshold to 
collapse. This observation may be useful to estimate the 
proximity to prompt gravitational collapse. Very near 
the threshold one thus may expect / 2-0 “S- /peak- In con¬ 
trast, /peak — /spiral typically increases with increasing 
Mtot, and above the threshold to collapse a spiral pat¬ 
tern during the dynamical collapse could still produce a 
weak peak in the GW spectrum, as in [s^ . 


IV. EMPIRICAL RELATIONS FOR DOMINANT 
AND SECONDARY PEAK FREQUENCIES 

For our sample of EoSs Fig. [5] shows /peak, /spiral 
and / 2-0 as a function of the compactness M/R of 
the nonspinning, individual NSs (at infinite separation) 
for Mtot = 2.7 Mq (with the compactness in units of 
c = G = I). We find strong correlations that can be 
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FIG. 7: /spiral VS. the compactness, but for different binary 
masses. Solid lines show empirical relations. The dashed line 
is taken from (see text for explanations). 



FIG. 8: / 2-0 vs. the compactness, but for different binary 
masses. Solid lines show empirical relations. The dashed line 
is taken from (see text for explanations). 


described by the following quadratic fits: 

/peak[kHz] = I99(M/i?)2 - 28.1(M/i?) + 2.33, (1) 
/spirai[kHz] = 358(M/i?)2 - 82.1(M/i?) + 6.16, (2) 
/ 2 -o[kHz] = 392(M/i?)2 - 88.3(M/i?)-h 5.95. (3) 

The maximum deviations of the data used for these fits 
are 140 Hz, 86 Hz and 153 Hz for /peak, /spiral and / 2 - 0 , 
respectively. If the compactness is determined from a 
measured frequency by inverting Eqs. these max¬ 

imum deviations imply errors of 3%, 3% and 4% in the 
compactess for /peak, /spiral and / 2 _o, respectively. (Note 
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that the slope in the relation for /spiral is somewhat flat¬ 
ter.) Thus, the accuracy of these empirical relations is 
similarly good for /peak and /spiral and only slightly worse 
for / 2 - 0 - Even if one assumes an error of 10% uncer¬ 
tainty in the frequency determinations from an actual 
GW observation, (which at least for /peak is too pes¬ 
simistic d^) this would add a 6%, 8% and 5% per cent 
error in the compactness determination. Similar empir¬ 
ical relations hold, with varying accuracy, for each total 
binary mass. As already pointed out and explained for 
/peak in [13 , even better empirical relations are obtained 
for 1.35-1.35 Mq binaries when the above three frequen¬ 
cies are examined as function of the radius of a nonro¬ 
tating NS with 1.6 Mq. Notice that all three empirical 
relations follow similar trends. 

Because the frequency of the /spiral peak agrees with 
the rotational frequency of the antipodal bulges near the 
surface of the remnant, a scaling of /spiral with the com¬ 
pactness is not unexpected, which explains the similari¬ 
ties with the behavior of /peak- In fact, for the models 
with Mtot = 2.7 Mq one finds very tight relations be¬ 
tween /spiral and /peak, and between / 2 _o and /peak- 

A measurement of /peak is already sufficient to accu¬ 
rately constrain the radius of nonrotating NSs (and thus 
the EoS) [H, m, 113 if the total mass has been obtained 
accurately from the inspiral signal, as it is likely to be 
the case for distances within which /peak has the re¬ 
quired signal-to-noise ratio to be detectable with second- 
generation interferometers (e.g. [13 113)- A detection of 
the weaker secondary peaks, with similar M/R depen¬ 
dence as /peak, could further optimize the constraints on 
the EoS. 

Figure [7] displays /spiral as a function of the compact¬ 
ness M/R oi the nonspinning, individual NSs for differ¬ 
ent EoSs and for different Mtot in the range 2.4-3.0 Mq 
(the most likely range of total binary masses [13)- The 
same figure, but for / 2 - 0 , is shown in Fig. [SI Especially 
for /spiral we discover that there exist tight relations be¬ 
tween the compactness and the secondary frequencies for 
fixed binary masses. 

In [H, [13 the existence of a single, universal, mass- 
independent relation (dashed line in Figs. [6] to [8]) between 
the frequency of the strongest secondary GW peak (de¬ 
noted there as /i) and M/R was proposed (there was no 
distinction of two different secondary peaks, as we find 
here). However, this result was based on using a limited 
set of five EoSs of soft or moderate stiffness (with cor¬ 
responding maximum masses of nonrotating NS only up 
to 2.2 Mq) as well as on different chosen mass ranges 
for each EoS with a spread of only 0.2 Mq in the total 
binary mass. 

In contrast to [H, [13, within our larger sample of 
ten EoSs (that includes stiff EoSs with maximum masses 
reaching up to 2.8 Mq) and for a more representative to¬ 
tal binary mass range of 2. 4-3.0 Mq (same for all EoSs), 
such a mass-independent, universal relation does not ex¬ 
ist (dashed curve in Figs. [5] to [5]). Figures |7| and |5|show 
that there is a large spread of up to 500 Hz both in the 


/spiral vs. M/i? relation and ill the /2-0 VS. M/i? relation 
for an expected range of total binary masses, and that 
the data cannot be described by a single function (dashed 
line). Even if one consistently chooses the strongest sec¬ 
ondary peak among /spiral and / 2-0 in each case, there 
does not exist a mass-independent universal relation with 
compactness. The absence of a mass-independent, uni¬ 
versal relation implies that the procedure for EoS con¬ 
straints and binary mass determinations as described 
in [13 [13 is not valid. Notice also that in Fig. [SI the 
relation proposed in [13, [13 describes either /spiral at 
low compactness or / 2-0 at high compactness, which is 
consistent with the expectations from our classification 
scheme describing which of the secondary peaks domi¬ 
nates the GW spectrum. However, for intermediate val¬ 
ues of the compactness the merger will be of Type H and 
both types of secondary peaks can be present with com¬ 
parable amplitude, which further complicates the defini¬ 
tion of a single /i frequency, as was assumed in [13 [33- 
We also point out that detections of the secondary peaks 
may possibly be less accurate since the peaks are broader 
in comparison to the main peak and sometimes do not 
stand out clearly. 

Instead of a universal mass-independent relation, we 
find that there exist useful empirical relations only for 
fixed binary masses (shown as thin line segments in 
Figs. 0 and j^), such as the case shown in Fig. [51 For 
merger events sufficiently close to allow a detection of 
post-merger GW peaks an accurate determination of the 
binary masses is expected (e.g. [13). With the up¬ 
coming GW detectors the dominant post-merger oscil¬ 
lation frequency has been shown to be detectable with 
high precision for mergers within several Mpc using a 
morphology-independent burst analysis [13. Significant 
improvements in the detectability are possible for more 
sophisticated search algorithms, e.g. matched filtering, 
which, however, require detailed modelling of the ex¬ 
pected signal features to which our present study should 
contribute by clarifying the relation between dominant 
and secondary GW peaks. Depending on the highly un¬ 
certain binary NS merger rate, a detection may succeed 
with Advanced LIGO/Virgo [13 or with discussed up¬ 
grades [43443 (the latter may yield a sensitivity increase 
at high frequencies of a factor three). Since the strength 
of the secondary features is somewhat lower than that 
of the dominant peak (in relation to the anticipated in¬ 
strument noise curves), direct detections of the secondary 
GW peaks can be expected for the planned Einstein Tele¬ 
scope [13 [13, EO, mi, 113 j unless the merger rate is on the 
more optimistic side as defined in [3, which may enable 
an earlier observation. The exact detection rate of sec¬ 
ondary peaks is hard to quantify not only because of the 
uncertain merger rate, but also because of the varying 
strength and prominence of the secondary peaks depend¬ 
ing on the exact model. 

In this study we focus on equal-mass binaries. To as¬ 
sess the impact of the binary mass ratio, we perform 
additional simualtions for 1.3-1.4 Mq binaries with the 


DD2 EoS [Hill, LS220 EoS H and NL3 EoS [Hill. 
The GW spectra in comparison to the symmetric case 
exhibit qualitatively the same features. The frequencies 
of the three different GW peaks are very similar to the 
equal-mass case with only small deviations of at most 2%, 
4% and 6% for /peak, /spiral and / 2 - 0 , respectively, com¬ 
pared to the results of the 1.35-1.35 Mq binaries. Also 
the evolution of the central lapse function shows qualita¬ 
tively the same behavior. This suggests that at least for 
moderately asymmetric binaries our classification scheme 
applies as well. It also implies that the frequency depen¬ 
dencies and their implications discussed above hold quan¬ 
titatively with a good accuracy, although the secondary 
peaks show somewhat larger deviations, which may im¬ 
pede their use for EoS constraints if the mass ratio is not 
measured accurately from the GW inspiral signal. 

V. SUMMARY AND OUTLOOK 

We present a unified picture of the post-merger dynam¬ 
ics and GW emission in binary NS mergers. Aside from 
the secondary GW peak produced by a nonlinear cou¬ 
pling between quadrupole and quasi-radial oscillations, 
we identify a new mechanism by which another secondary 
peak in the GW spectrum appears: it is caused by a spiral 
deformation, the pattern of which rotates slower with re¬ 
spect to the double-core structure of the inner remnant. 
Based on the presence and strength of these two sec¬ 
ondary peaks, we introduce a classification scheme of the 
post-merger dynamics and of the GW spectrum. More¬ 
over, the new mechanism allows for an explanation of 
low-frequency modulations seen in a number of physical 
characteristics of the remnant, such as the central lapse 
function, the maximum density, the separation between 
the two cores and the GW amplitude. 

We find that for fixed total binary masses the frequen¬ 
cies of the secondary GW peaks show a tight relation with 
the compactness of nonrotating, individual NSs. These 
relations follow similar trends as the relation between the 
dominant peak frequency and the NS compactness (or, 
equivalently, the NS radius, as shown in l^)- For 


the dominant peak and for the secondary peak associ¬ 
ated with the orbital motion of the antipodal bulges we 
find the relations between frequency and compactness to 
be similarly tight. We rule out the existence of a uni¬ 
versal, mass-independent relation for secondary peaks as 
proposed in [H, . 

Identifying the type of a merger event and especially 
the determination of the quasi-radial frequency by /o = 
/peak — / 2-0 can lead to an estimate of the threshold 
mass for black-hole collapse and thus to an estimate of 
the maximum mass of nonrotating NSs (see [l^). The 
insights from our classification scheme and from the fre¬ 
quency dependencies found here can provide a basis for 
constructing future GW template waveforms, increasing 
the chan ges for the observation of the post-merger GW 
emission (40|. 

We will further investigate whether the frequency dif¬ 
ferences between the dominant and secondary peaks can 
clarify the nature of a detected secondary peak or if de¬ 
tailed comparisons between spectra of different binary 
masses are needed. We will also explore unequal-mass 
binaries (anticipating a strong impact of the mass-ratio 
on /spiral) and analyze the relevance of our classification 
scheme for the mass ejection and torus formation of NS 
mergers and for accompanying phenomena, such as r- 
process nucleosynthesis [59|, l60|, electromagnetic coun¬ 
terparts (bll - l^ and short gamma-ray bursts [13, [13■ 
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